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ABSTRACT 

We present a model of inflation based on the anomaly-induced effective action of gravity in 
the presence of a conformally invariant Hilbert-Einstein term. Our approach is based on 
the conformal representation of the fields action and on the integration of the corresponding 
conformal anomaly. In contradistinction to the original Starobinsky's model, inflation can be 
stable at the beginning and unstable at the end. The instability is caused by a slowing down of 
inflation due to quantum effects associated to the massive fermion helds. In supersymmetric 
theories this mechanism can be linked to the breaking of SUSY and suggests a natural way 
to achieve graceful exit from the inflationary to the FLRW phase. 



Introduction 

Inflation Q automatically solves five of the six basic cosmological problems : 1) the monopole 
problem, 2) the horizon problem, 3) the flatness-curvature-entropy problem, 4) the rotation problem, 
and 5) the large-scale homogeneity versus small-scale inhomogeneity problem. The minimum number 
of e-folds of inflation required can be (roughly) estimated in many different ways. As an example, 
take a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) model, then the scale factor a = a(t) in the 
matter-dominated (MD) and radiation-dominated (RD) eras evolve as a ~ t 2 / 3, and a ~ t x l 2 respectively. 
Since at present to ~ I0 18 sec (15 Gy) and a(to) ~ 1.5 x 10 10 lyr ~ 10 28 cm, it follows that the scale 
factor at the end of the RD epoch (t ~ 10 12 sec ~ 10 5 yr) was an = a(to)/zR = 10 24 cm, where 
zr ~ a(to)/aji = (l0 18 /10 12 ) 2 ^ 3 = 10 4 is the redshift at that time. From this the scale factor at the 
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Planck time (~ 10 — 44 sec) should be a* P = (l0~ 44 /10 12 ) an ~ 10 cm, which is of course untenable! 
Therefore, to make this number to match up the correct Planck length, ap ~ I0~ 33 cm, we need to 
supplement the standard FLRW evolution with an early inflation period in which the number of e-folds 
of inflation should be around £ ~ 65: 

e « = ^ = 10 2 9^£~67. (1) 
ap 

At present, however, the classical cosmological problems that motivated inflation are no longer regarded as 
the strongest motivation for inflationary cosmology. For example, the relation "homogeneity — ► flatness" 
is not true. Nowadays we have homogeneous open models of inflation. Besides providing the solution to 

1 Talk given at SUSY 2002: the 10th International Conference on Supersymmetry and Unification of Fundamental 
Interactions, DESY, Hamburg, Germany, June 17-23; also at GRG 11: International Conference on Theoretical and 
Experimental Problems of General Relativity and Gravitation, Tomsk, Russia, July 1-8 2002. 



some cosmological problems, the important thing at present is that inflationary models can be (and will 
be) more and more accurately tested (something that one could not suspect 10 years ago) through their 
specific predictions on the metric and density perturbations, which should be consistent with structure 
formation and the anisotropies of the CMB H . These are nowadays the real facts behind £ > 65 in 
Eq.(|l]), rather than the previous and similar heuristic argumentations. 

Inflation, however, does not solve the sixth cosmological problem, the cosmological constant (CC) 
problem |Q. The fact that the CC has been measured non-zero and positive |5| poses a new challenge and 
may require novel approaches. One possibility is to think of the renormalization group (RG) evolution 
of the cosmological parameters ||. We will see that this same approach can be applied to the study of 
cosmological inflation, if we identify the RG scale with the expansion rate p = H(t). This identification 
can be understood as follows. At low energy the dynamics of gravity is defined by Einstein's equations 

R^u - - Rg^u = SttG n (T m „ + 5 M „ A ) , (2) 

where Gn — 1/Mp is Newton's constant. Let us use the value of the curvature scalar (R) to construct 
an order parameter for the gravitational energy. By dimensional analysis the RG scale p for gravity 
is naturally associated with i? 1 / 2 . From Eq.(|J) we see that this is equivalent to take p ~ y/Tfi /Mp. 
But in the cosmological setting the basic dynamical equations refer to the scale factor a(t) of the FLRW 
metric, and so we must re-express the graviton energy in terms of it. The 00 component of @ yields the 
well-known Friedmann-Lemaitre equation 

o f a\ 2 8n . , . k , . 

H =(-) =TT^ ^ + A 3 
\a ) 3 Alp a z 

The space curvature term can be safely set to zero (k = 0). The spatial components of (p|), combined 
with the 00 component (|J), yields the following dynamical equation for a(t): 

j(p + 3p- 2 A) a. (4) 
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In these equations p = pu + PR is the total energy density of matter and radiation, and p is the pressure. 
In the modern Universe p ~ and p ~ p° M . Moreover, from the recent supernovae data (5j, we know that 
A and p° M have the same order of magnitude as the critical density p° c . Therefore, the source term on 
the r.h.s. of (||) is characterized by a single dimensional parameter yf p®/M'p, which according to Eq. (||) 
is nothing but the experimentally measurable Hubble's constant Hq. This is obviously consistent with 
the expected result sjT^jM'p in the general case because Tff ~ p° M ~ p° c for the present-day universe. 
Therefore, we conclude that the Ansatz 

fi - R 1 / 2 - H(t) (5) 

is reasonable and we assume that this identification takes place at each stage of the cosmological evolution. 
With this guiding principle, a semi-classical description of gravitational phenomena of quantum matter 
in a curved classical background should be possible. In particular, in the early universe the value of H (t) 
decides which matter particles are active degrees of freedom for the RG evolution of the parameters. 
Therefore, particles whose masses satisfy M > H(t) will be decoupled from the relevant quantum effects 
at time t. Clearly, if one would be able to concoct a natural mechanism by which H(t) progressively slows 
down after the universe has achieved a sufficient number of e-folds of inflation (namely £ > 65), then 
inflation should eventually stop and the FLRW regime could perhaps start. While many authors have 
looked for a suitable scalar field (so-called "inflaton") capable of realizing this scenario Q, in the original 
Starobinsky's model inflation was attempted by looking for a self-consistent solution of Einstein's 
equations when they are modified to include the vacuum quantum effects. Unfortunately, in Starobinsky's 
model inflation is unstable from the very beginning (with the flat space stable), so that one has to fine- 
tune the initial conditions to insure £ > 65 before inflation stops B . It would be much desirable to find an 
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improved framework where the self-consistent solution appears first as a stable inflationary solution (hence 
independent of the initial conditions, even though space-time is unstable) and such that subsequently 
(after £ > 65 is fulfilled) inflation becomes unstable and the universe transits into the stable and flat FLRW 
space-time. Indeed, a modified Starobinsky's model like that is possible, provided that we can arrange the 
condition H(t) — > and at the same time distort the stability regime thanks to a change in the number 
of active degrees of freedom, e.g. due to a phase transition from a supersymmetric Grand Unified Theory 
(GUT) into the Standard Model (SM) of the strong and electroweak interactions. Such a scenario has been 
first discussed in Ref. || . It is based on a modification of the anomaly-induced effective action |?], |8[ [l0| 
resulting from a prior full conformization of the classical action for gravitational fields [jllj , including the 
Hilbert-Einstein term, and matter fields [jl2[ , and on the decoupling of the supersymmetric particles at 
low energy |l3[. Furthermore, there are strong indications that the spectrum and the amplitude of the 



gravitational waves in this model |t4|, 15 are in agreement with the existing CMBR data|3|]. Recently, 
also the stability with respect to small perturbations of the conformal factor of the metric has been 
studied in the presence of a cosmological constant || [H| . 

2. Local conformal invariance and effective action 

The expansion of an homogeneous, isotropic universe means a conformal transformation of the metric 
<Jnv(t) ~ * ° 2 (??) <?ai!; j where a(r]) = exp a(rf) and r] is the conformal time (drj = dt/a(rj)). Suppose that 
one starts from the conformal invariant formulation of the Standard Model of the strong and electroweak 
interactions |T^ , |l6[ and of gravityjll], and then one uses the well-known methods to derive the anomaly- 
induced action p"8| . It is natural to think that the latter can be applied at high energies, where 
the masses of the matter fields are negligible. At the classical level, the theory which results from 
this procedure is always equivalent to the original theory. Nevertheless, in the quantum theory the 
equivalence is destroyed by the anomaly, which can be calculated explicitly. In particular, the massive 
fields may also contribute to it. Besides the anomalous terms, there are the conformal invariant quantum 
corrections to the classical vacuum action. Our first purpose is to construct such a formulation of the SM 
in curved space-time which possesses local conformal invariance in d = 4. Actually, the procedure can be 
applied to any gauge theory, e.g. the SM and extensions thereof, including GUT's and supersymmetric 



generalizations like the Minimal Supersymmetric Standard Model (MSSM) 19 . 

The original action of the theory includes kinetic terms for spinor and gauge boson fields, as well 
as interaction terms, all of them already conformal invariant. As for scalars (e.g. Higgs bosons) we 
suppose that their kinetic terms appear in the combination g^d^ipd^ip + 1/6 • Rip 2 providing the local 
conformal invariance. The non-invariant terms are the massive ones for the scalar and spinor fields, but 
also the Hilbert-Einstein term giving General Relativity at low energies. In all these cases the conformal 
non-invariance is caused by the presence of dimensional parameters m 2 H , m, Mp = l/G. The central 
idea is to replace these parameters by functions of some new auxiliary scalar field x- F° r instance, we 
replace |fT^] 

2 m H 2 ™ 2 Mp 2 

m H^ W X, M P ^— X , (6) 

where M is some dimensional parameter, e.g. related to a high scale of spontaneous breaking of dilata- 
tion symmetry |jl2f . Then the scalar and fermion mass terms become quartic interactions and Yukawa 
couplings respectively, 

lJd 4 xV^^V 2 x\ f*xJ=g%H, X , (7) 

which are of course (local and global) conformal invariant. Furthermore, the Hilbert-Einstein term gets 
conformized too: 

S*eh = ~ 167r g M 2 / d4x ^9 [Rx 2 + & {d X f } . (8) 
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After setting \ — > M this expression becomes identical to the ordinary gravitational term, and from 
the ordinary mass terms for scalars and fermions are recovered at the same time . This fixing can 
be called "conformal unitary gauge" in analogy with the unitary gauge of ordinary gauge theories, and 
the scale M can be associated to the vacuum expectation value of the spontaneously broken dilatation 
symmetry at high energies fl2]l 

It is supposed that the new scalar field \ takes the values close to M, especially at low energies. 
But, there is a great difference between \ and M with respect to the conformal transformation. The 
mass does not transform, while \ does. Then, the action of the new model becomes invariant under the 
conformal transformation 



(9) 

which is performed together with the usual transformations for the other fields 

g^^g^e 2 *, ip^ V e-\ ^->^e- 3 / 2CT . (10) 

Thus, in the matter sector our program of "conformization" is complete. When we quantize the theory, 
it is important to separate the quantum fields from the ones which represent a classical background. In 
order to maintain the correspondence with the usual formulation of the SM, we avoid the quantization 
of the field x which will be considered, along with the metric, as an external classical background for the 
quantum matter fields. It is well known (see, e.g. |18| ) that the renormalizability of the quantum field 
theory in external fields requires some extra terms in the classical action of the theory. The list of such 
terms includes the nonminimal term of the J Rip 2 -type in the Higgs sector, and the action of external 
fields with the proper dimension and symmetries. The higher derivative part of the vacuum action has 
the form 

S vac = J d A x^g~ {hC 2 + l 2 E + l 3 V 2 R } , (11) 

where, £1,2,3 are some parameters, C 2 is the square of the Weyl tensor and E is the integrand of the 
Gauss-Bonnet topological invariant. Now, since there is an extra field \, the vacuum action should be 
supplemented by the ^"dependent term. The only possible, conformal and diffeomorphism invariant, 
terms with dimension 4 are (Q) and the J x 4 -term. The last contributes to the cosmological constant, 
which we suppose to cancel and do not consider here. Its effect is reported elsewhere || [l5| . 

The next step is to derive the conformal anomaly in the theory with two background fields and 
X • The anomaly results from the renormalization of the vacuum action including the terms @ and ( |ll|) . 
For the sake of generality, let us suppose that there is also some background gauge field with strength 
tensor F^ v . Then the conformal anomaly has the form 

<T; >= - {wC 2 + bE + cV 2 R + dF 2 + f [R X 2 + 6 (d X ) 2 }} , (12) 

where w, b, c are /3-functions for the parameters l\, l 2 , h', f is the /3-function for the dimensionless 
parameter l/(16irG M 2 ) of the action (||), and d is the /3-function for the gauge coupling constant. The 
values of w, b and c depend on the matter content (JVj being the number of particles with spin i): 

_ N + 6N 1/2 + 12N 1 N + llN 1/2 + 62N 1 _ N + 6N 1/2 - 18N r 

W = 120 • (4tt) 2 ' = 360 • (4tt) 2 ' ° = 180 • (4tt) 2 ' (13) 

Recall that the condition for stable inflation is c > M. Then one can play with various models. For 
instance, from the previous equation it follows that the particle content of the SM (No — 4,Ni/ 2 — 
24, Ni = 12) leads to c < (unstable inflation), which suggests that with only SM matter fields the 
inflation period cannot be easily sustained, and it might be insufficient 0. On the other hand, for the 
MSSM [[HI (No = 104, N 1/2 = 32, Ni = 12) one has c> (stable inflation) etc. Clearly, we need physics 
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beyond the SM in order to possibly arrange a graceful transit from a regime of stability (insuring the 
condition £ > 65) into another of instability that might hopefully end into the FLRW phase. 

We will argue that the presence of the non-zero /3-function / could be the necessary quantum dynam- 
ical mechanism for the graceful exit. From direct calculation using the Schwinger-DeWitt method (see 
e.g. @) we get 

J ^ 3 (4tt) 2 M 2 ' 1 ' 

where N% are the number of Dirac spinors with masses mi. We note that bosons do not contribute to /. 

In order to obtain the anomaly-induced effective action, we put g^ u = g^ v ■ e 2a and X = X ' e ~ a ■> 
where the metric g^ has fixed determinant and the field x does not change under the conformal 
transformation. Then, the solution of the equation for the effective action T proceeds in the usual way 
. Disregarding the conformal invariant term we arrive at the following expression || : 

f = J d i x^j{wC 2 + b(E - ^V 2 i?) + 2baA + dF 2 + 

+ /[i?X 2 + 6(9x) 2 ]}--^^ Jd^x^-gR 2 . (15) 



3. The role of masses in slowing down inflation 

In order to understand the role of the particle masses in the anomaly-induced inflation, let us consider 
the total action with quantum corrections 

St = Smatter + SeH + S vac + T . (16) 

One of the approximations we made was to disregard higher loop and non-perturbative effects in the 
vacuum sector. Another approximation is that we take only the leading-log corrections. Usually, this is 
justified if the process goes at high energy scale. If the quantum theory has UV asymptotic freedom, 
the higher loops effects are suppressed, and our approximation is reliable. At the low-energy limit, we 
suppose that the massive fields decouple and their contributions are not important. Then Eq. ( |l6| ) can 
be presented in the form 

a = / rf4 ^ { ( - + f° ) [ Rx 2 + e m 2 ] 

— ^ — — G?er ^ F 2 j + S ma tter + high, deriv . terms . (17) 

One can see that the modifications with respect to the case of free massless fields |Q are an additional 
/-term and the contribution to anomaly due to the background gauge fields. In order to restore the 
Hilbert-Einstein term and get the inflationary solution, we fix the conformal unitary gauge and put 
X = X e " — M. Furthermore, wc CcLn choose the conformably flat metric Q^iu — Vf^v- 

Then the gravitational 

part of the action ( |l7|) becomes 



Sgrav — 
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J d 4 x { 2b (d 2 a) 2 - (3c + 26) [(da) 2 + d 2 a)} 2 - 

^-Tr'Hi-*)**}- < 18 > 
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Figure 1 . (a) Plot of a — In (a) versus the physical time t as a result of the numerical analysis 
of Eq.|Tg[); t is given in units of 16 ir/Mp and we fixed the parameter ( p?C[ ) as / = 1CP 4 . Initial 
data: a(0) = 1 , a(0) = Hi , a(0) = H 2 , 'a(O) = Hf . In this time interval, inflation does not 
stop, yet; (b) As in (a), hut extending the numerical analysis until reaching an approximate 
plateau marking the end of stable inflation. 



Computing the equation of motion of a = In a in terms of the physical time t (where dt = a(rj)drf) we 
find 

2 .... ., .... /'_ 1'A .2.. ..2 Ml ( .2\ 

a a + iaaa — ( 5 H a a + aa — a a + aa 

C I 07TC V 



2/Af 2 



In a 
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a a + aa 



2/Af 2 



dF 2 
6ca 



0. 



(19) 



An exact solution of ( |Ts| ) does not look possible, but it can be easily analyzed within the approximation 
that / is not too large. Then the new terms (collected in the second line of Eq. (|H^)) can be considered 
as perturbations. Moreover, the last two of them are irrelevant, because during inflation they decrease 
exponentially with respect to the other terms. Thus, in this approximation, the only one relevant change 



is the replacement Mp - 
dimensionless parameter 



Mp 



l-/Ino(t) 



where for future convenience we have introduced the 



~ Mp ^ 3tt Af 2 



(20) 



Notice that / is given by Eq. ([II|) and so / does not depend on the scale M. Since / is small, the effect of 
the masses may be approximated through the modification of the inflation law a(t) = ao e Hlt according 

ton 



#i = 



Mz 



x/-16tt6 



1/2 

l-fhxa(t) =H(t), 



(21) 



To substantiate our claim, we have solved Eq. ( |19| ) directly using the numerical methods. The plots 
corresponding to the numerical solution of the Eq. are shown in Fig. 1. Since in the first period of 
inflation masses do not play much role and the stabilization of the exponential inflation proceeds very 
fast, the initial data (in both Eq. ( |2l| ) and the plots of Fig. 1) were chosen according to the exponential 
inflation law. According to the numerical analysis, the total number of e-folds in the "fast phase" of 
inflation (until the Hubble constant becomes comparable to the transition scale M* where instability 
develops) is about 10 for our particular values of the parameters, and at the last stage the expansion 
essentially slows down. 



2 We remind the reader that the coefficient b is negative for any particle content, see 
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Figure 2. Plot of H(t) — a{t)/a{t) versus t as a result of the numerical analysis of Fq.(pl|) 
and parameter values as in Fig.l: (a) H(t) near the onset of the plateau; (b) H(t) well over 
the plateau. 

The chosen value of the parameter (|2^) / = 1CP 4 in the plot is, as we warned before, independent of 
the scale M, and it determines where the process of stable inflation finishes as well as the number of 
e-folds. On the other hand, if the scale M* is chosen near the typical SUSY GUT value Msusy ~ 
Mx ~ 10 16 GeV for supersymmetry breaking at high energies, some of the spinor masses m, will be 
of order Msusy ~ 10 16 GeV. The latter assumption is indeed sound because the mi will include the 
supersymmetric fermions associated to the super-heavy gauge and Higgs bosons at the GUT scale, and 
so mi ~ Mx ~ Msusy- From Eq. (|2(]) it follows that the parameter / will be numerically smaller 
than the one we have assumed in Fig. 1, and consequently the amount of inflation will be larger. But 
the important qualitative point is that for any value of / the approximate plateau eventually appears 
and signals the end of stable inflation. Also notice from Fig. 1 that the initial evolution is close to the 
exponential inflation, but after that the expansion slows down due to the quantum effects of massive 
fermions. 

4. Graceful exit from anomaly-induced inflation 

Recall from Eq.(|^) that H(t) sets the scale of the RG running for the gravitational part. So if we 
consider the SUSY breaking and the corresponding change in the number of active degrees of freedom, 
then the necessary and sufficient condition for the applicability of our approach is that H (t) decreases 
from the initial value aboutQ Mp/^-Wirb ~ 10 18 GeV, down to the lower scale H = M* < M SUSY . 
The outcome is that the evolution according to (HJ) lasts until reaching the scale M* , and after that most 
of the SUSY particles are decoupled, the inflationary solution becomes unstable and the FLRW phase 
can start. In fact, the crucial point is the existence of a nonvanishing / as it eventually tempers stable 
inflation allowing favorable conditions for the universe to tilt into the FLRW phase. 

To justify our claim that M* < Msusy, recall that for a really successful exit from the inflationary 
phase we need to make sure that the amplitude of the gravitational waves is consistent with the observable 
range of anisotropy in the CMBR. This will be the case if during the last 65 e-folds of the inflation, 
the expansion rate H(t) does not exceed 10~ 5 Mp. Then the fluctuations in the amplitude h of these 
waves, 5h/h = H/Mp, will preserve the measured fluctuations in the temperature of the relic radiation 
according to the relation 5h/h = 5T/T — 0(1O~ 5 ). At the lowest end of the inflation interval this 
condition corresponds, in our framework, to fix the instability value M* « 10~ 5 Mp = 10 14 GeU. It 
means that, in reality, we expect that after the onset of the approximate plateau in Fig. I (b), where 
SUSY breaking occurs, the universe will take a while before entering the FLRW phase, i.e. the latter will 
actually initiate at some point well over the plateau where H = A/* ~ 10~ 5 A/p. To better assess this 
issue we have numerically analyzed H(t) over the plateau, see Fig. 2. We see that H(t) decreases very 

3 Notice that | 16tt6| = 0(1) in the MSSM, and it is much larger than 1 in any typical SUSY GUT. 
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fast on it. For instance, from the comparison of Fig. 2(a) and Fig. 2(b), we find that a 15% increase of 
the time after the onset of the plateau amounts H (t) to diminish two orders of magnitude. So in general 
H{t) will decrease further below Msusy, and the difference between M* and Msusy at the moment of 
the transition can be significant, say one or two orders of magnitude. Hence Msusy can be 10 16 GeV and 
this does not create problems with the CMBR. Let us also notice the oscillatory behavior of H(t) after 
reaching the plateau, i.e. when the system is about jumping into the FLRW phase. Interestingly enough, 
this behavior could perhaps be related to the reheating phenomenon, which is of course indispensable 
before the universe stabilizes in the FLRW regime. 

Overall, we arrive at a consistent picture of the graceful exit in this inflationary scenario. No inflaton 
field is needed, but only the dynamical work of gravity itself, provided one starts from a renormalizable, 
fully conformal-invariant classical picture, together with the trace anomaly of the matter fields at the 
quantum level. Moreover, according to (fhf), the obtained picture is universal, for it does not depend on 
the choice of the dilatation symmetry breaking scale M. If interpreted physically, one can put constraints 
on M using the macroscopic forces mediated by the field a, demanding that this forces should have the 
sub-millimeter range, similarly as in ]2(i| ]. 
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